Draft version February 1, 2008 

Preprint typeset using LATp^X style emulateapj v. 04/03/99 



MEASURING ANGULAR DIAMETER DISTANCES THROUGH HALO CLUSTERING 

ASANTHA COORAY 12 , WAYNE Hu 1 , DRAGAN HUTERER 3 , AND MlCHAEL JOFFRE 1 

1 Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637 
3 Department of Physics, University of Chicago, Chicago, IL 60637 
E-mail: (asantc,whu,dhuterer,joffre)@oddjob. uchicago.edu 

Draft version February 1, 2008 
ABSTRACT 

Current and upcoming wide-field surveys for weak gravitational lensing and the Sunyaev-Zel'dovich 
effect will generate mass-selected catalogues of dark matter halos with internal or followup photometric 
redshift information. Using the shape of the linear power spectrum as a standard ruler that is calibrated 
by CMB measurements, we find that a survey of 4000 deg. 2 and a mass threshold of 10 14 Mq can be used 
to determine the comoving angular diameter distance as a function of redshift. In principle, this test also 
allows an absolute calibration of the distance scale and measurement of the Hubble constant. This test is 
largely insensitive to the details of halo mass measurements, mass function, and halo bias. Determination 
of these quantities would further allow a measurement of the linear growth rate of fluctuations. 

Subject headings: cosmology: theory — large scale structure of universe 



1. INTRODUCTION 

A number of observational efforts are now underway or 
being planned to image the large-scale structure of the 
universe spanning a range of redshifts. These wide- field 
surveys typically cover tens to thousands of square degrees 
of the sky: the ongoing Sloan Digital Sky Survey (SDSS), 
the weak gravitational lensing shear observations with in- 
struments such as SNAP or the Large Aperture Synoptic 
Survey Tel escope, and surveys o f the Sunyaev-Zel'dovich 
effect (SZ; |Sunyaev k Zel'dovich| 1980). 

In addition to their primary science goals, these surveys 
are expected to produce catalogues of dark matter halos, 
which in the case of lensing a nd SZ surveys ar e expe cted to 



be ess entially mass selected (Wittman ct al. 2001; Holder 



et alj 2000). Lensing and other optical surveys are par 



ticuiarly promising in that they will provide photo metric 
redsh ifts on the member galaxies of a given halo (e.g. , Hogg 
et al. 1998); this will render accurate determination of the 
halo redshift. Halo number counts as a function of redshift 
is a well-known and powerful cosmological test. Here we 
consider the additional information supplied by the angu- 
lar clustering of halos. 

A feature in an angular power spectrum of known phys- 
ical scale and originating from a known redshift can be 
used to measure the angular diameter distance between us 
and this redshift; this has most notably been applied to 
the case of the cosmic microwave background (CMB) to 
determine the distance to redshift z ~ 10 3 . The angular 
power spectrum of halos provides a similar test based on 
the standard ruler defined by its shape. In the adiabatic 
cold dark matter model for structure formation, this stan- 
dard ruler is essentially the horizon at matter-radiation 
equality and its absolute physical scale can be directly cal- 
ibrated with CMB anisotropy data. In principle then, one 
can determine the angular diameter distance as a function 
of redshift and test the properties of the dark energy. 

As a purely geometric test, this method is largely in- 
sensitive to uncertainties in the the mass function and the 
relationship between the halo masses and the actual ob- 
servables, e.g. the SZ temperature decrement or lensing 
aperture mass. The bias of the halos is scale-dependent 




Fig. 1. — Angular power spectrum of halos with M > 10 Mq in 
a wide-field survey in bins of z = 0.3 — 0.4 and 1.0 — 1.2. The binned 
errors are 1-<t, and assume a survey of 4000 deg 2 , within reach of 
upcoming weak lensing and SZ surveys. The angular power spec- 
trum at high redshifts is shifted towards the right proportional to 
the increase in the comoving angular diameter distance. The oscil- 
lations are due to baryons but we ignore the additional information 
they contain. 



only on small (non-linear) scales, and can in principle be 
extracted to arbitrary precision from A-body simulations. 
If and when these quantities are securely known, one can 
extract further information from the amplitude and small- 
scale behavior of the power spectrum. In particular, the 
linear growth rate and non-linear scale provide extra han- 
dles on the dark energy. 

For illustrative purposes, we adopt the ACDM cosmol- 
ogy with energy densities (relative to critical) of Q, m — 0.35 
in matter, f2f, = 0.05 in baryons, f^A = 0.65 in vacuum 
energy, the dimensionless Hubble constant of h — 0.65, 
and a scale-invariant spectrum of primordial fluctuations, 
normalize d to the present day galaxy cluster abundance 
(a 6 = 0.9; |vTana fc Liddle] 1999). 
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2. ANGULAR POWER SPECTRUM 

The angular po wer spec trum of halos in ith redshift bin 
is a Limber (1954; Kaiser 1992) projection of the the halo 
number density power spectrum 



r;= i dzW?(z)^Lp, 



hh 



1 



(1) 



where Wi(z) is the distribution of halos in a given redshift 
bin normalized so that J dzWi(z) = 1, H(z) is the Hub- 
ble parameter, and (La is the angular diameter distance 
in comoving coordinates. Note that Wi{z) comes directly 
from the observations of the number counts as a function 
of redshift and depends on the mass function and mass 
sensitivity of the employed observable. 
If the halos trace the linear density field, 



P hh (k;z) = (b M r(z)D'(z)P lm (k;0), 



(2) 



where (6m) is the mass-averaged halo bias parameter, 
P hu (k; 0) is the present day matter power spectrum com- 
puted in linear theory, and D(z) is the linear growth func- 
tion S lin (k;z) = D(z)S l[n (k; 0). A scale-independent halo 
bias is com monly assumed in the so-called "halo model" 
(e.g. Seljak 2000) and should be valid at least in the linear 
regime. Equation (|l]) then becomes 

Cj = J dzWf{z)F{z)P^iL ] Q \ , (3) 
H(z) 
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(4) 



The underlying linear power spectrum contains two 
physical scales: the horizon at matter radiation equality 



k cq = \/2fl m ifo(l + z cq) oc Q m h 



(5) 



which controls the overall shape of the power spectrum, 
and the sound horizon at the end of the Compton drag 
epoch, k s (il m h 2 , fli,h 2 ), which controls the small wiggles in 
the power spectrum. The angular or multipole locations 



of these features shift in redshift as L 



We propose the following test: measure C\ in several red- 
shift bins and, using the fact that l eq scales with dj^{zi), 
constrain the angular diameter distance as a function of 
redshift. To be conservative, we ignore the additional in- 
formation suppled by l s . 

In Fig. |l|, we illustrate the proposed test. The two 
curves show the halo power spectra in two redshift bins: 
0.3 < z < 0.4 and 1.0 < z < 1.2. The angular power spec- 
trum corresponding to the higher redshift bin is shifted to 
the right in accordance to the ratio of angular diameter 
distances (81/1 ~ 53 a/ (1a)- Much of this shift simply re- 
flects the Hubble law, dA ~ z/Hq. Since the physical scale 
of the two features — the overall shape of the spectrum 
and the baryon oscillations — can be calibrated from the 
morphology of the CMB peaks, these measurements can 
in principle be used to determine the Hubble constant in- 
dependently of the distance ladder and distance to last 
scattering surface. 

In addition to the horizontal shift due to the change 
in angular diameter distance, the power spectra in Fig. Py 



are shifted vertically due to the change in F(z) (Eq. B). 
By ignoring the information contained in F(z), this pure y 
geometric test is robust against uncertainties in the mass 
selection, mass function and linear bias. Of course, if these 
uncertainties are pinned down independently, both F(z) 
and the halo abundance in Wi(z) will help measure the 
growth rate of structure. 

3. PARAMETER ESTIMATION 

Even though the angular diameter distance test is ro- 
bust against uncertainties in the halo selection function, 
number density and bias, these quantities enter into the 
consideration of the signal-to-noise for a realistic survey. 
We will focus on a survey of 4000 deg. 2 with a detection 
threshold in mass of 10 14 Mq out to z = 2. We use a to- 
tal of 9 bins in redshift; note that the cluster photometric 
redshift accuracy is expected to be much smaller than the 
bin width. The mass threshold is consistent with those ex- 
pected from upcom ing le nsing and SZ effect surveys (see 
Kruse fe Schneider] 1999; |Holder et aT] 2000; Joffre et al. 
in preparation) and the survey area is consistent with a 
planned SZ survey from the South Pole Telescope (Carl- 
strom, private communication). To compute Wf(z) we 
adop t the predictions of the Press-Schechter mass function 
(PS; Press fe Schechtei| 1974). This mass funct ion, along 
with the halo bias prescription of Mo fc White (1996), is 
also used to predict the mass-averaged halo bias (6m) (z)- 

Assuming Gaussian statistics, we can express the uncer- 
tainty in the measurements of the angular power spectrum 
as 



AC} = 



{Cf + Nf 



VG + l/2)/sky ' 



(6) 



where / s k y = 9^^/129600 is the fraction of the sky cov- 
ered by a survey of dimension 0d C g in degrees and Nf is 
the noise power spectrum. We assume that the dominant 
source of noise is the shot-noise so that Nf = 1/iVj, where 
Ni is the surface density of the halos in the zth redshift 
bin. We use the PS mass function to predict Ni. In Fig. [j], 
the two bins contain roughly ~ 4 and 6 halos/deg. 2 above 
our minimum mass. In the same Figure, we show band 
power measurement errors following Eq. |6j. 

To estimate how well halo clustering can recover cosmo- 
logical information, we construct the Fisher matrix 
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(l + l/2)/ Bky dCf dCj 



(7) 



where a and (3 label parameters that underly the power 
spectra. Since the variance of an unbiased estimator of 
a parameter p a cannot be less than (F _1 ) QCe , the Fisher 
matrix quantifies the best statistical errors on parameters 
possible with a given data set. 

We choose l mm = 27r/0do g when evaluating equation (|^) 
as it corresponds roughly to the survey size. The precise 
value does not matter for parameter estimation due to the 
increase in sample variance on the survey scale. Given 
our crude Gaussian approximation of the shot-noise, we 
choose a conservative i^ ax corresponding to the multipole 
at which the noise and sample variances are equal Nf — 
C? : 'max ranges from 200 at low redshift bins to 400 at 
high redshift. At low redshifts this cutoff is slightly in the 
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Fig. 2. — (a) The errors (1-cr) on angular diameter distance as 
a function of redshift. We have binned the halos in 8 redshift bins 
between and 1.6. The larger errors are with no prior assumption 
on the cosmological parameters that define the transfer function 
while the smaller errors are with MAP (Temp) and Planck (Pol) 
priors. In (b), we show relative errors in the distance. 



non-linear regime, but at redshifts greater than 0.8 or so 
one is well within the linear regime. Therefore, such a 
low l max largely eliminates uncertainties in the modeling 
of scale-dependent halo bias in the non-linear regime. 

Due to the dependence of C\ on P lm (fc), all cosmo- 
logical parameters that change the shape of the matter 
power spectrum across the scales probed by halos also af- 
fect the measurement of distance. The shape of the trans- 
fer function is determined by fl m h 2 and fl^h 2 while the 
overall slope is determined by a scalar tilt n s . Note that 
these parameters will be accurately determined from CMB 
anisotropy observations. Since the CMB peaks probe the 
same range in spatial scale as the halo power spectrum, 
one is relatively insensitive to deviations from a pure ini- 
tial power law. 

When estimating expected errors on distance, we will 
consider several sets of priors on these cosmological pa- 
rameters. These priors follow Table 2 of Eisenstein et al 



(1999) and correspond to constraints expected from MAP 
and Planck with and without polarization. Though baryon 
oscillations contain cosmological information, in order to 
be conservative against possible non-linearities in the bias 
we ignore the information present in the bar yon oscilla- 
tions, and employ the smooth fitting function of Eiscnstcir] 
& Ht|(1999). Our results are then very weakly dependent 
on the fiducial value or priors on £li,h 2 . In case baryonic 
features in the angular power spectrum are detected, we 
expect additional cosmological information to be gained 
using the proposed test. 

4. RESULTS & DISCUSSION 

We first consider the measurement of the angular diam- 
eter distance d z A = d,A(Zi). In addition to the cosmological 
parameters fl m h 2 , fl^h 2 and n s , we include a set of pa- 
rameters F % = F(zi) which allow the normalization of the 
i power spectra to float independently. Both d\ and F l 



approximate underlying functions as piecewise flat across 
each bin. 

In Fig. ||(a), we show three sets of errors: the largest er- 
rors assume no prior knowledge on the transfer function, 
while the smaller errors correspond t o priors fr om MAP 
(Temp) and Planck (Pol) respectively (Hu et al. 1999). In 
Fig. ||(b), we show the fractional percentage errors on the 
distance. The errors in the lowest bins tell us how well 
one can estimate the Hubble constant, while the slope of 
dji(z) around z ~ 1 provides information on cosmology. 
dA is best determined around z « 0.7 at the level of 10%. 
For systematic errors in the redshifts to dominate the error 
budget, the inferred mean redshift of the bin must differ 
from the true value by ~ 10%, which is far above the er- 
rors we expect for individual halos even from photometric 
techniques at z < 1. Precision at higher redshift is not 
required because of the large statistical errors. The errors 
on distance estimates d\ are correlated at the 5% level due 
to remaining uncertainties in parameters that affect all d l A 
(e.g., £l m h 2 ). 

In a realistic cosmology dA(z) is smoothly varying. Since 
fl m h 2 is already taken as a parameter, we parameterize 
dA with the Hubble constant h = _ffo/100km s _1 Mpc^ 1 
and the equation of state of the dark energy w, the ratio of 
pressure to density, assuming a flat Universe. We bin halos 
following the binning scheme in Fig. ||. In Fig. ^a we show 
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Fig. 3. — The errors (all 1-cr) on (a) h and ui and (b) Q m and w, 
using distance information only. In both (a) and (b) we show errors 
for halos with priors following MAP (Temp) and Planck (Pol). In 
(b), for comparison, we also show errors on D.™ and iij from CMB 



(Planck with temperature and polarization (Hu et al 
type la SNe with SNAP mission. 



1999)) and 
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that a strong degeneracy in h and w remains even with 
Planck priors because (1a is only accurately recovered in a 
small redshift range. Of course, an external determination 
of h would break this degeneracy. 

An alternate way of breaking the degeneracy is to em- 
ploy other cosmological probes of Q m and w. As shown 
in Fig. [j](b), different linear combinations of fl m and w 
will be determined by halos and the CMB due to the dif- 
ferent redshift ranges probed. Although each constraint 
alone may not be able to pin down w, halos and the CMB 
combined allow very interesting constraints even under our 
conservative assumptions. 

Uncertainties in the mass threshold and the scale de- 
pendence of halo bias are potential caveats to these con- 
clusions. A mass threshold that differs from the assumed 
10 14 Mq value would not bias the angular diameter dis- 
tance results since they only utilize redshift and power 
spectrum shape information. However it would affect the 
errors due to a rapid decrease in the number density of 
halos with threshold mass: at 4 x 10 14 Mq the error on w 
increases increase by a factor of ~ 3. 

A scale-dependent bias that can be predicted actually 
aids in the determination of angular diameter distances: 
the scale-dependence acts as another standardizable ruler 
for the test. Indeed, the scale-dependence of the bias as a 
function of halo mass is something th at can be precisely 
determined from A-body simulations ( Kravstov fc KlypTq 
1999). A more subtle problem is introduced by the ad- 
dition of uncertainties in the mass threshold or selection 
function. Since the bias is also mass-dependent, the uncer- 
tainty in the mass threshold SM/M — 0.1 translates into 
the uncertainty in the mass-averaged bias S (6m) / (6m) = 
0.03. To investigate a scale dependent bias, we model it 
as 



b M (k, z) = (6 M ) (z) 



1 + / 



P nl (fc;z) 



- 1 



(8) 



P lin (fc; z) 

where / is a dimensionless parameter meant to interpo 




Fig. 4. — The 1-a error on w as a function of the prior on bias. 
The four curves assume Planck (Pol) priors on Q m h 2 , n s and Q b h 2 
to define the linear power spectrum. The solid line is with no prior 
on In A, and h. The dotted line includes a prior of 0.2 in In A, while 
the long dashed line is with an addition prior of 0.1 in h. The dot- 
dashed line is a highly optimistic scenario with exact A, a prior of 
0.1 in h, and using halo angular power spectrum information out to 
l max of 1000, deeply within the non-linear regime. 



late bias between the linear (/ — > 0) and the non-linear 
(/ — > 1) regimes. Note that in the halo approach to clus- 
tering, the non-linear mass power spectrum is a sum of the 
halo power spectrum and contributions due to dark matter 
within halos; therefore, the halo power spectrum cannot 
be larger than the non-linear power spectrum. Taking a 
fiducial model with / = and adopting MAP (Temp) pri- 
ors, we find that marginalizing over / increases the error 
on w by less than ten percent. 

To the extent that the mass threshold, halo bias, and 
mass function are known, the amplitude of the halo power 
spectra can be used to measure the linear growth rate. We 
conclude by estimating the level at which these quantities 
must be determined to yield additional constraints on w. 
In Fig. we plot the marginalized errors on w as a func- 
tion of the assumed fractional prior on (6 M ) for various 
independent constraints on A and h. Since cosmological 
information captured in linear growth is determined by 
relative amplitude variations in C\, the knowledge of the 
overall normalization A is not crucial. For example, go- 
ing from no prior knowledge of A to a Gaussian prior with 
width of 20% of the fiducial value of A results in a decrease 
in cr(w) of ~ 25%. Errors on the mass selection function 
bias the measure of w. Using an extension to the Fisher 
matrix approach, we determined that a 25% systematic 
offset in mass threshold from the fiducial value of 10 14 
Mq leads to systematic bias in w of 0.05 from its fi ducial 



value of -1.0. Lcnsing si mulations (Metzler et al. 1999 



Reblinsky & Bartlemann 1999) indicate that the calibra- 
tion of projection effects at this level will be challenging 
but feasible to achieve. 

Clearly future surveys which can identify dark matter 
halos as a function of redshift contain valuable information 
beyond the evolution of their number abundance. As the 
theoretical modeling of the halo distribution and empirical 
modeling of the selection process improve, the correlation 
function of the halos can provide not only the angular 
diameter distance, but also direct measurements of the 
growth of large-scale structure. 
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